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Abstract. This paper discusses the overall quality of lUE 
high resolution data processed with the NEWSIPS soft- 
ware in terms of flux and wavelength accuracy. It also de- 
scribes the processing of NEWSIPS high resolution spec- 
tra within the framework of the ESA "lUE Newly Ex- 
tracted Spectra" (INES) System. This system provides the 
lUE high resolution data in two formats: the high resolu- 
tion "concatenated" spectra, in which the spectral orders 
are connected, eliminating the overlap regions through a 
procedure designed to optimize the signal-to-noise ratio at 
the edges of the orders, and the "rebinned" spectra, i.e. 
the high resolution concatenated spectra resampled into 
the low resolution wavelength domain. Our study reveals 
the existence of a significant discrepancy in the wavelength 
scales of short and long wavelength NEWSIPS high reso- 
lution spectra. The INES processing applies a correction 
of -t-17.7 km s^^ to the wavelength scale of the high reso- 
lution SWP spectra in order to provide an internally con- 
sistent velocity scale. Similarly, suitable corrections have 
been applied to long wavelength small aperture spectra. 
The new wavelength scale is, within the errors, in good 
agreement with the optical velocity scale. 

Key words: Methods:Data Analysis; Techniques: Image 
Processing; Ultraviolet: General 



1. Introduction 

The lUE NEWSIPS piocessing system was developed with 
the aim of creating a "Final Archive" of lUE data to 
be made available to the astronomical community as a 
legacy after the end of the project. This archive would 
include all the lUE spectra re-processed with improved 
algorithms and up-to-date calibrations. The NEWSIPS 
processing system is fully described by Nichols and Linsky 



Send offprint requests to: R. Gonzalez-Riestra 
* Previously: ESA-IUE Observatory 
** Affiliated to the Astrophysics Division, SSD, ESTEC 
Correspondence to: ch@laeff.esa.es 



(1996) and Nichols (1998). Technical details are given in 
the NEWSIPS Manual (Garhart et al. 1997). The intro- 
duction of new techniques to perform the geometric and 
photometric corrections led to a substantial improvement 
in the signal-to-noise ratio of the final spectra. Further im- 
provements in the quality of the high resolution data arise 
from the new method to determine the image background 
and from the improved ripple correction and absolute cal- 
ibration (Cassatella et al. 1999, hereinafter Paper II). The 
background subtraction has always been one of the more 
critical issues in the processing of lUE high resolution 
spectra, particularly at the shortest wavelengths, where 
orders crowd and an accurate estimate of the background 
is essential for a correct flux extraction. This problem has 
been overcome in NEWSIPS through the derivation of a 
bi-dimensional background (see Smith 1999 for a descrip- 
tion of the method) . 

The goal of the INES processing system was to cor- 
rect the deficiencies found during the scientific evalua- 
tion of the data processed with NEWSIPS for the lUE 
Final Archive, and to provide the output data to the 
users in a simple way requiring a minimum knowledge 
of the operational and instrumental characteristics of lUE 
(Wamsteker et al. 1999). The modifications introduced in 
the processing of low resolution data have been described 
by Rodriguez-Pascual et al. (1999, Paper I). As for high 
resolution data, the INES system provides two output 
products derived from the NEWSIPS MXHl (i.e. high res- 
olution extracted spectra) files: the "concatenated" spec- 
trum, where the spectral orders are merged eliminating 
the overlap regions, and the "rebinned" spectrum, which 
is the concatenated spectrum resampled to the low reso- 
lution wavelength step. Both concatenated and rebinned 
spectra include an error vector calibrated in absolute flux 
units. The inconsistency between the high resolution short 
and long wavelength scales in NEWSIPS has been cor- 
rected for in the INES concatenated spectra. 

In the flrst part of this paper we discuss the over- 
all quality of NEWSIPS high resolution spectra in terms 
of accuracy, stability and repeatability of wavelength and 
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Table 1. Radial velocities obtained from NEWS IPS high resolution spectra 



Target 



SWF 



LWP 



LWR 



RRTel 
COph 

BD+28 4211 
HD 60753 
HD 93521 
BD+75 325 
A Lep 
C Gas 
ri UMa 



-69.5±6.5 [106] em. lin. 
-24.7±3.9 [60] SII, Sill 
-22.8±6.4 [130] SII,SiII,CII 
18.8±6.6 [62] SII,OI,SiII 
-38.8±4.0 [68] SII,OI,SiII 
-16.4±4.6 [121] SII,SiII,CII 



-49.3±3.0 [170] em. lin. 
-13.4±2.6 [24] FeII,MnII 
-5.3±3.6 [40] Mgll K 
32.4±6.9 [70] Mgll K 
-20.7±6.9 [6] Mgll K 
6.2±2.2 [43] Mgl 
21.2±3.9 [37] Mgll K, Mgl 
-1.1±4.2 [9] Mgll K 
-1.1±3.1 [17] Mgll K 



-51.0±4.4 [132] em. lin. 
-10.2±2.5 [69] Fell, Mnll 
-1.1±4.8 [21] Mgll K and H, Mgl 
29.5±5.8 [24] Mgll K and H, Mgl 
-19.9±5.4 [18] Mgll K and H 
6.6±3.1 [33] Mgll K 
19.5±4.9 [44] Mgll K and H 
1.5±4.4 [72] Mgll K and H 
3.1±7.1 [82] Mgll K and H 



The entries for each star are: mean radial velocity, rms (km s^^), number of measurements, and lines used (em. lin.: emission 

lines) . 

All measurements refer to the large aperture. 



flux measurements (Sectional). The second part deals with 
the INES processing of high resolution data, describing 
the order concatenation and rebinning procedures (Sec- 
tions 3.1 and 3.2, respectively). Finally, the application of 



the correction to the wavelength scale is discussed (Section 



2. NEWSIPS Data quality evaluation 

The overall quality of lUE high resolution spectra pro- 
cessed with NEWSIPS has been evaluated by studying 
the accuracy and the stability of wavelength determina- 
tions, the accuracy of equivalent width measurements, the 
flux repeatability and the residual camera non-linearities. 
The analysis is based on a large number of spectra, 
mainly of the lUE standard stars. The spectra have been 
corrected for the echelle blaze function and calibrated in 
terms of absolute fluxes according to the procedure de- 
scribed in Paper IL Hereinafter, wavelengths are assumed 
to be in the heliocentric reference frame and in vacuum. 

2.1. Wavelength accuracy 

To assess the wavelength accuracy three aspects have been 
considered separately: 

a) the accuracy and repeatability of wavelength mea- 
surements of a given spectral feature in several spectra of 
the same star, 

b) the stability of the wavelength scale along the full 
spectral range, 

c) the consistency of radial velocity determinations ob- 
tained from the SWP, LWP and LWR cameras. 

2.1.1. Expected accuracy 

One of the most important limitations to the wavelength 
accuracy of lUE high resolution spectra are the target ac- 
quisition errors at the nominal center of the spectrographs 



entrance apertures. These, if large enough, can also affect 
the quality of the ripple correction (see Paper II). 

From the NEWSIPS dispersion constants and the 
central wavelengths of the spectral orders it can be 
readily deduced that the velocity dispersion corre- 
sponding to one pixel on the image is practically constant 
all through the range covered by the cameras, and namely: 

AV=7.73 ± 0.05 km s"! for SWP 
AV=7.26 ± 0.09 km s^^ for LWP 
AV==7.26 ± 0.03 km s^^ for LWR 

Taking into account that the plate scales are 1.530, 
1.564 and 1.553 arcsec/pix for SWP, LWP and LWR, re- 
spectively (Garhart et al. 1997), an acquisition error of 1 
arcsec along the high resolution dispersion direction would 
lead to a constant velocity offset of 5.1 km s^^ for SWP, 
4.6 km s-i for LWP and 4.7 km s^^ for LWR. Since the 
pointing/tracking accuracy is usually better than 1 arc- 
sec, we can consider 5 km s^^ as a reasonable upper limit 
to the expected wavelength accuracy. Wavelength errors 
substantially larger might arise internally in the data ex- 
traction procedures. 

2.1.2. Repeatability of wavelength determinations 

To obtain a reliable information on the self-consistency of 
wavelength determinations, we have attempted to reduce 
the effects of spectral noise by performing a large num- 
ber of measurements of selected narrow and symmetric 
absorption lines from the interstellar medium, which are 
strong in some of the lUE calibration stars, as well as of 
many emission lines in RR Tel. Out of the lUE standards, 
we have selected BD-l-28 4211, HD 60753, HD 93521, 
BD-f75 325, A Lep (HD 34816), C Gas (HD 3360) and rj 
UMa (HD 120315). In addition, we have also used spectra 
of the star ( Oph (HD 149757). The present measurements 
refer only to large aperture spectra. 
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The interstellar lines selected for the SWP range were: 
SII 1259.520 A, Sill 1260.412 A, 01 1302.168 A, Sill 
1304.372 A and CII 1334.532 A, 1335.703 A. For the long 
wavelength range we have used the K and H components of 
the Mgll doublet at 2796.325 A and 2803.530 A and Mgl 
2852.965 A. In the case of ^ Oph we have also measured 
Mnll 2576.877 A, 2594.507 A and several Fell lines. The 
cases in which the distortion of the profile due to close-by ' 
reseau marks (in particular for 01 1302 and CII 1334) pre- 
cluded the accurate determination of the line position have 
not been taken into account in the final statistics. Labo- 
ratory wavelengths have been taken form Morton (1991). 

The mean values of the radial velocities, the corre- 
sponding rms deviation and the number of independent 
measurements are reported in Table |^ for each target and 
camera. According to this table, the rms repeatability er- 
ror on radial velocities, averaged over the three cameras 
is 4.6±1.5 km s^^. This value is smaller than the upper 
limit of about 5 km s~^ expected from acquisition errors. 
Considering the presence of spectral noise, we can safely 
conclude that the repeatability of wavelength (or velocity) 
determinations is satisfactory. 

The results in Table ^ indicate that the radial veloc- 
ities derived from the two long wavelength cameras are 
consistent while, on the contrary, the velocities derived 
from SWP spectra are systematically more negative. This, 
and other considerations about the consistency of radial 
velocity determinations from the three cameras will be 
discussed in Section ^. 

2.1.3. Stability of the wavelength scale along the full 
spectral range 

We have studied the accuracy of the wavelength scale over 
a wide spectral range to look for possible time-dependent 
distortions across the camera faceplate. To this purpose, 
we have selected 6 SWP, 11 LWP and 3 LWR spectra 
of the emission line object RR Tel obtained at different 
epochs. For each spectrum we have measured the peak 
wavelengths of several emission lines chosen among those 
reasonably well exposed and with the cleanest profiles, 
covering the full spectral range. The highest excitation 
lines, such as those from [MgV], were purposely excluded 
because they provided systematically higher negative ra- 
dial velocities probably due to stratification effects within 
the nebular region. The mean radial velocities of RR Tel 
are -69.5±6.5 km s'^, (SWP), -49.3±3.0 km s'^ (LWP), 
and -51.0±4.1 km s^^ (LWR). The total number of mea- 
surements are 106, 170 and 132 for SWP, LWP and LWR, 
respectively. Since the errors are of the same order as the 
repeatability errors quoted in the previous section, we con- 
clude that the wavelength scales do no present appreciable 
distortions over the wavelength range covered and, within 
the observational errors, are stable over the period of time 
considered (1983-1994, 1985-1995 and 1978-1983 for SWP, 
LWP and LWR, respectively). 
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Fig. 1. A portion of the longest available SWP exposure 
of RR Tel (SWP20246: 820 minutes). The two Lyman a 
lines arise from geocoronal emission filling the large and 
the small apertures. The emission features labeled with an 
asterisk are due to overspilling of the strong NV emission 
into the adjacent order. 



2.1.4. Radial velocity determinations from the Mg II 
doublet 

The present analysis has revealed the existence of an in- 
consistency in the radial velocities derived from the Mgll 
K (2796.32 A) and H (2803.53 A) lines as measured in 
the LWP camera, where the two lines are present in both 
orders 82 and 83. To quantify this discrepancy, we have 
measured the velocities of the Mg II interstellar lines in 
89 spectra of five lUE standard stars. 

We find that, in the LWP camera, the radial velocity 
difference Vh-Vk is -10.8±2.5 km s"-^ when measured in 
order 83, and 11.7±1.5 km s~^ when measured in order 82. 
We find also that there is a discrepancy in the velocity of 
the K line measured in the two orders (V83-V82=-20.5±0.9 
km s~^), while the velocities of the H line measured in 
the two orders are consistent within 2 km s^^, on average. 
Since the velocity derived from other interstellar lines (e.g. 
Mgl 2852.965 A) is fully consistent with the measurements 
from the K line in order 83, we conclude that only this line 
provides correct radial velo city values. In INES concate- 
nated spectra (Section |3.l|) the K line comes from order 
83 and the H line from order 82. Therefore, there is a sys- 
tematic difference Vk-Vh=-8.8±1.3 km s~^ between the 
velocities determined from the two lines, being the correct 
value only that given by the K line. 

A similar study performed in LWR spectra, where the 
K line appears only in order 83, shows that this problem 
is not present in this camera, where the two Mg II lines 
provide consistent velocity values: VK(order 83)-VH(order 
82)=-1.0±1.2kms^i. 
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2.2. Background extraction 

It has been repeatedly pointed out that the background 
extraction for high resolution spectra processed with 
lUESIPS was not accurate enough especially shortward 
of 1400 A in the SWP camera and 2400 A in the long 
wavelength cameras, as denoted by the negative fluxes as- 
signed to the wings of the strongest emission lines and 
to the core of the saturated absorption lines. As shown 
below, this effect is not present anymore in spectra pro- 
cessed with NEWS IPS, which makes use of an upgraded 
background determination procedure (Smith 1999). 

Overestimating or underestimating the background 
level leads to underestimating or overestimating the fluxes 
of the emission lines and the equivalent widths of the ab- 
sorption lines. In the following, we report the tests done on 
the accuracy of the equivalent widths to verify the correct- 
ness of the background extraction. In addition, we have 
used the repeatability of the equivalent widths determina- 
tions as an indirect test of the stability of the background 
levels. 




0.1 0.15 

E.W. Copernicus 



2.2.1. SWP 

Figure || shows the profiles of the NV doublet emission 
and of the broad Lyman a feature in the longest expo- 
sure available of RR Tel. It is clearly seen that the wings 
of these lines are not assigned negative values. Particu- 
larly interesting are the NV "ghost" lines marked with an 
asterisk, which are still present in NEWSIPS data, but 
sensibly fainter than in the lUESIPS spectra, most likely 
due to the optimized extraction slit. The presence of such 
spurious lines has recently been reported by Zuccolo et 
al. (1997) and ascribed to overspilling of the strong NV 
doublet into adjacent orders. 

To verify the accuracy of the background subtraction, 
we have compared the equivalent widths of the strongest 
interstellar lines in four spectra of C, Oph with those re- 
ported by Morton (1975) obtained from Copernicus data. 
These latter determinations are presumably not affected 
by background determination problems, unlike the lUE 
echelle spectra near the short wavelength end of the cam- 
eras. The results of the comparison are given in Fig. g and 
in Table y , which provides the mean and the standard de- 
viation of the four measurements. As it appears clearly 
from the figure, the NEWSIPS measurements are consis- 
tent, within the errors, with the values from Copernicus 
data, suggesting that the background evaluation for SWP 
spectra is essentially correct. 

The stability of the background subtraction has been 
evaluated by measuring the equivalent width of several 
strong interstellar lines in a large sample of spectra of two 
standard stars. The repeatability of the equivalent widths 
ranges from 10% for the strongest lines to 30% for the 
faintest ones (Table ^) . 



Fig. 2. Comparison of the equivalent widths of interstellar 
lines in the spectrum of ^ Oph as measured in Copernicus 
and in NEWSIPS spectra. Copernicus measurements have 
been taken from Morton (1975). The dashed line in each 
panel indicates the 1:1 relation. 

Table 2. Equivalent widths (in mA) in SWP spectra of (^ 
Oph compared with Copernicus measurements 



Line 


NEWSIPS 


Morton (1975) 


SIT 1250.59 


81±12 


100±3 


SII 1253.81 


92±11 


106±2 


SII 1259.52 


93±18 


112±1 


Sill 1260.42 


167±26 


170 


CI 1277.24 


62±13 


74 


01 1302.17 


221±20 


201±7 


Sill 1304.37 


115±13 


132±1 


CI 1328.83 


49±7 


52±6 


CII 1334.53 


193±17 


189±2 


CII 1335.70 


147±21 


140±4 


Sill 1526.71 


170±15 


190: 


Fell 1608.46 


120±15 


85: 


Sill 1808.01 


54±8 


87±11 


Aim 1854.72 


49±5 


57: 


Aim 1862.79 


28±2 


34: 



means uncertain value 



2.2.2. LWP 



We tested the accuracy of the background extraction from 
the cores of strongly saturated absorption lines and the 
wings of strongly saturated emission lines. Fig. shows 
a portion of a spectrum of SN1987A centered around the 
Mg II doublet at 2800 A. It appears from the figure that 
the cores of the absorption lines do not become system- 
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Table 3. Equivalent widths (in mA) in LWP and LWR 

spectra of C. Oph compared with Copernicus measure- 
ments 



Line 


LWP 


LWR 


Morton (1975) 


2.0x10-1° 




Znll 2026.16 


111±26 


117±32 


100: 




Fell 2382.76 


251±6 


246±11 


238±20 




ifU 


Mnll 2576.88 


133±12 


124±20 


138±21 


1 l.SxlO-lO 




Fell 2586.65 


225±10 


200±8 


214±19 






Mnll 2594.51 


112±8 


125±10 


116±5 


l.OxlO-lO 


- 


Fell 2600.17 


251±5 


253±6 


226±16 






Mnll 2606.48 


104±14 


120±8 


111±6 


5.0xlO~" 


- 




: means uncertain value 

Table 4. Repeatability of equivalent widths measure- 
ments 

SWF 



Line 



BD-h28 4211 BD-f 75 325 



SII 1250.59 32±9 (29) 34±8 (26) 

SII 1259.52 48±10 (28) 70±12 (35) 

SII 1260.42 92±12 (37) 128±14 (25) 

CII 1334.53 115±14 (36) 125±19 (35) 



LWP 


Line 


BD+28 4211 


BD-l-75 325 


Mgll 2796.53 
Mgl 2852.97 


257±20 (37) 
276±40 (11) 


369±50 (33) 
45±16 (10) 


LWR 


Line 


rj UMa 


C Gas 


Mgll 2796.35A 
Mgll 2803.53A 


65±11 (41) 
75±23 (41) 


354±22 (36) 
345±23 (36) 



Equivalent widths in mA. The number of measurements is 
given in brackets. 



atically negative, as expected for a correct background 
extraction. Shown in the same figure is another exam- 
ple, that of the strongly saturated Mg II emission dou- 
blet in the longest LWP exposure available of RR Tel 
(LWP25954): the lines wings do not reach negative fluxes, 
as required. 

A second test has been performed in four spectra of 
C Oph, in which we have measured the equivalent widths 
of some strong interstellar lines, and compared them with 
the Copernicus values given by Morton (1975). The results 
are summarized in Fig. and Table 0, where are given the 
mean and the standard deviation of the four spectra. 

Finally, we have measured the equivalent widths of 
the Mgll and Mgl lines in a large sample of spectra of 
BD-F28 4211 and BD-h75 325. The repeatability errors 
range from 35% for the faint Mgl line in BD-l-75 325, to 
less than than 10% for the strong lines. Results arc shown 
in Table |. 




v U^ 
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Fig. 3. Top panel: The region of the Mgll doublet in 
the spectrum of SN 1987A (LWP10194). Bottom panel: 
The same region in the deepest LWP exposure of RR Tel 
(LWP25954). 



2.2.3. LWR 

As for the LWR camera, we have verified the accuracy of 
the background subtraction by measuring the equivalent 
widths of six strong Fell and Mnll interstellar lines in four 
spectra of C Oph, and compared them with measurements 
based on Copernicus data. As shown in Fig. g, there is a 
good agreement between the two sets of equivalent widths, 
and no systematic departures are found. 

The repeatability of equivalent width determinations 
has been determined measuring the equivalent widths of 
the Mgll doublet in a large sample of spectra of rj UMa 
and ^ Cas and A Lep. The repeatability error ranges from 
30% for the weak lines to less than 5% for the strongest 
ones (Table |). 

2.3. Flux repeatability 

Two different tests have been performed to assess the flux 
repeatability of high resolution spectra. In the first one 
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Table 5. Flux repeatability of high resolution spectra 2.3.1. SWP 



SWP 



Band 


Order 


set A 


set B 


1185-1190 


116 


3.2 


3.5 


1225-1230 


112 


2.3 


4.3 


1285-1290 


107 


2.6 


3.2 


1375-1380 


100 


1.8 


3.0 


1485-1490 


93 


1.6 


3.9 


1595-1600 


86 


3.5 


3.5 


1795-1800 


77 


1.6 


2.5 



Set A: BD+28 4211: a) 4 spectra (Jun 86-Dec 86); b) 6 spec- 
tra (Nov 90- Jul 91), BD+75 325: 8 spectra (Dec 88-May 89), 
HD 60753: a) 4 spectra (Apr 80-May 80); b) 8 spectra (Feb 
88-May 89), HD 93521: 11 spectra (March 1987) 
Set B: BD-f 28 4211: 45 spectra (Dec 82- Aug 95), BD+75 425: 
53 spectra (Sep 85- Feb 95), HD 60753: 38 spectra (Jul 79- 
Mar 95) 

LWP 



Band 


Order 


set A 


set B 


2117-2122 


109 


4.9 


6.3 


2457-2462 


94 


2.0 


3.5 


2922-2927 


79 


1.8 


2.7 


3117-3122 


74 


3.0 


4.0 



9 



Set A: BD-h75 325: a) 11 spectra (March 85-Dec 85); b) 

spectra (Jan 89-Dec 89); c) 8 spectra (Feb 92-Nov 92) 

Set B: BD-h28 4211: 34 spectra (Dec 82-Aug 95), BD-h75 325: 

37 spectra (Sep 85-Feb 95), HD 60753: 24 spectra (Jul 86- Aug 

95) 

LWR 



Band 


Order 


Set A 


Set B 


2117-2122 


109 


4.2 


5.5 


2457-2462 


94 


3.4 


6.3 


2922-2927 


79 


1.8 


4.3 


3077-3082 


75 


2.4 


5.0 



Set A: ( Gas: a) 4 spectra (Dec 81-Dec 82); b) 5 spectra (Sep 
83- Jul 84), ri UMa: a) 6 spectra (Mar 81- Apr 82); b) 5 spectra 
(Aug 82- Jun 86) 

Set B: C Cas: 30 spectra (Feb 81-Feb 87), r] UMa: 42 spectra 
(Sep 78- Jul 90) 



we selected restricted samples of spectra obtained close 
in time, and measured the ripple corrected net fluxes, i.e. 
without applying the time sensitivity degradation and the 
temperature dependence corrections. The second test has 
been performed on larger samples covering extended pe- 
riods of time, measuring the absolutely calibrated fluxes, 
which include time and temperature corrections. In all 
case we have averaged the flux over a narrow wavelength 
interval free of lines. The flux repeatibility was defined 
as the percent rms deviation from the mean value. The 
results are summarized in Table ^. 



For the first test we have used 41 high resolution spec- 
tra of lUE calibration standards grouped into sets of data 
with a similar exposure level and obtained close enough in 
time. The test was done in six bands 5 A wide. In Table || 
(under "set A" ) we report the percent rms deviation. The 
repeatability of spectra obtained sufficiently close in time 
is about 2%. The second test was made on a larger number 
of spectra with similar exposure time, without restricting 
the date of observation ("set B"). This test provides re- 
peatability errors ranging from 3 to 4%. These somewhat 
larger errors are due to the intrinsic uncertainties of the 
sensitivity degradation correction algorithm. 



2.3.2. LWP 

The tests performed are similar to those described above 
for the SWP camera. The spectra in "Set A" consist of 
three groups each containing images obtained in a re- 
stricted period of time. The flux repeatability was eval- 
uated in four wavelength bands 5 A wide. In this cam- 
eras, the repeatability errors can reach the 5% level near 
the short wavelength end of the camera, but are a fac- 
tor of two lower in the central bands. A similar test per- 
formed on a larger set of spectra needing correction for 
the time-dependent sensitivity degradation ( "Set B" ) pro- 
vides errors slightly larger, conflrming that the sensitivity 
degradation algorithm adopted for the LWP camera is es- 
sentially correct. 



2.3.3. LWR 

The flux repeatability was evaluated in four wavelength 
bands 5 A wide. The test performed on spectra taken close 
in time, "Set A" , shows that the repeatability errors reach 
4% near the short wavelength end of the camera, decreas- 
ing in the region of maximum sensitivity and increasing 
again at the longest wavelengths. In the "Set B" spectra 
the repeatability is worse, reflecting the uncertainties in 
the time degradation correction and also the instability of 
the camera after it ceased to be routinely used. 



2.4- Flux Linearity 

Despite the linearity correction applied during the pro- 
cessing, residual non linearities are still present in lUE 
data. This effect has been evaluated in Paper I for low 
resolution data. In what follows we discuss this effect in 
high resolution spectra. The method followed consists on 
studying a set of spectra of the same star with different 
exposure times obtained, whenever possible, very close in 
time (preferably on the same observing shift) and with 
similar camera temperatures. The variation of the flux 
with the level of exposure (or the exposure time) defines 
the flux linearity. Unfortunately there exist few sets of 
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Table 6. Linearity study for high resolution spectra 



SWP 



-6% at 1185 A to +4% at 1785 A for the 50% exposure 
and up to -5% for the 70% exposure. 



Image t/t(opt) 1185A 1285A 1485A 1785A 2 4 2 LWP 



41467 
41346 
42309 
41435 
41466 
42260 
41495 



0.50 
0.68 
0.70 
0.90 
1.00 
1.00 
1.80 



0.94 
0.92 
0.97 
0.96 
0.98 
1.02 
0.96 



0.95 
0.98 
0.93 
0.99 
0.98 
1.02 
0.98 



0.97 
0.95 
0.98 
1.01 
1.00 
1.01 
0.98 



1.05 
1.02 
1.00 
0.99 
0.98 
1.03 
0.99 



LWP 



Group t/t(opt) 2 120 A 2460A 2925A 3132A 



a 
b 
c 
d 

6 
f 



0.27 
0.41 
0.67 
0.83 
1.00 
1.33 
2.07 



1.03 
1.00 
0.99 
0.95 
1.00 
0.93 
1.01 



1.03 
1.02 
1.02 
0.98 
1.00 
1.03 
1.04 



0.98 
0.99 
1.03 
1.00 
1.00 
1.00 
1.04 



0.98 
0.98 
1.00 
0.97 
1.00 
1.03 
1.08 



a: 2 spectra of BD+28 4211 

b: 1 spectrum of BD+28 4211 

c: 2 spectra of BD+28 4211 

d: 2 spectra of BD+75 325 

e: reference group: 4 spectra of BD+75 325 and 4 spectra of 

BD+28 4211 

f: 1 spectrum of BD+75 325 

g: 1 spectrum of BD+75 325 



LWR 



Image 


t/t(opt) 


2120A 


2460A 


2925A 


3080A 


9955 


0.50 


0.75 


0.98 


0.99 


0.91 


9954 


1.00 


1.00 


1.00 


1.00 


1.00 


9113,9953 


2.25 


1.18 


1.05 


1.00 


1.01 


8116 


2.50 


1.12 


1.04 


1.01 


1.00 



high resolution data suitable for this study, and a slightly 
different approach has been taken here. The results are 
summarized in Table ^. In general, these results are in 
good agreement with those derived in Paper I. 



There is not any complete set of high resolution data of 
the same star which allows to study the LWP camera lin- 
earity. We have constructed average spectra of different 
exposure levels of the two standard stars BD+28 4211 
and BD+75 325, and divided them by the correspond- 
ing 100% spectrum. The exposure levels covered range 
from 27% to 207% of the optimum exposure time. The 
test was performed in four wavelength bands 5 A wide, 
selected for being relatively free from strong absorption 
lines. The maximum departures from linearity, reaching 
8%, are found for the 133% level at 2120 A, and for the 
207% level near the regions of maximum sensitivity of the 
camera. The latter deviation can be easily understood in 
terms of saturation. 

2.4.3. LWR 

The LWR high resolution linearity test has been per- 
formed with five images of the standard star HD 93521 
obtained in the period July 1980 to February 1981, cover- 
ing the range of exposure times from 50% to 250% of the 
optimum value. The test was performed in four bands 5 
A wide. The maximum departures form linearity (up to 
25%) are found at the shortest wavelengths, where fluxes 
are underestimated by 25% for the 50% exposure and over- 
estimated by 12% for the 250% exposure. 



3. INES Processing of High Dispersion Spectra 

The starting point for the INES processing of high resolu- 
tion data are the NEWSIPS MXHI files. The spectra have 
not been re-extracted from the bi-dimensional files, as in 
the case of the low resolution data (see Paper I). The 
INES system provides two output spectra for each high 
dispersion image: the "concatenated" and the "rebinned" 
spectra. Both include a modified wavelength scale and the 
error vector calibrated in absolute flux units. All these as- 
pects are discussed in the following sections. 



2.4.1. SWP 

The most complete set of SWP spectra appropriate to 
assess the linearity consists of seven images of the white 
dwarf CD-38 10980 obtained in the period April-August 
1991, with exposure times ranging from 50% to 180% of 
the optimum value (200 min). For these images we have 
measured the mean flux in five 5 A wide bands. The fluxes 
in each band have been averaged together and divided by 
the mean flux in the two 100% exposures. The ratios so 
obtained indicate departures from linearity ranging from 



3.1. The high resolution concatenated spectra 

The main features of the INES concatenated spectra are: 

a) The overlap regions between adjacent orders are 
suppressed in such a way that the less noisy portion of 
the orders is retained. 

b) The error vector is calibrated in absolute flux units. 

c) The wavelength scale is modified to make consis- 
tent the radial velocities obtained from the short and the 
long wavelength cameras. The wavelength sampling of the 
MXHI files has been retained. 
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Fig. 4. Example of the concatenation procedure for the 
SWP camera. The top panel shows the individual orders 
from the MXHI file. The bottom panel shows the final 
concatenated spectrum and the calibrated errors (dotted 
line) . 



d) The spectra are provided as FITS tables having 
the same format as the low resolution spectra, i.e. they 
contain only four columns: wavelength, absolute flux, er- 
ror and quality factor. This setting reduces significantly 
the download time for remote data retrieval, and simpli- 
fies considerably the structure of the NEWS IPS MXHI 
files, which contain additional information, not relevant 
for most investigations, such as the position of the orders 
on the bi-dimensional image, the height of the extraction 
slit and the order number for each extracted point. 

The procedure followed to concatenate the spectral or- 
ders is described in detail in the next paragraphs. 

3.1.1. The concatenation procedure 

The most critical point in the concatenation of adjacent 
echelle orders is a suitable definition of the "cut wave- 
lengths" so that only the highest quality data points of 
the overlap region are retained. 



Fig. 5. As figure Q for the LWP camera. 

In lUE high resolution spectra, the signal-to-noise 
level at the edge of the orders is different in the short 
and long wavelength cameras. In the LWP and LWR 
cameras the signal-to-noise is always lower at the short 
wavelength end of the orders (except in the highest orders, 
see below). The contrary happens in the SWP camera, 
where the long wavelength edge of the orders is much 
noisier than the short wavelength one. Taking this into 
account, the cut wavelengths have been defined as follows: 

SWP: 

^cut — ^start + \^end ~ ^start)/-^- (1) 

LWP/LWR: 

AcMt = ^start + 2 X [Xend ~ Xstart)/^- (2) 

where 

Xcut'- cut wavelength between orders m and m-1 
Xstart'- start wavelength of overlap region (order m-1) 
Xend '■ end wavelength of overlap region (order m) 

These expressions are valid for all spectral orders ex- 
cept order 125 in the LWP camera and orders 120 to 125 
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Table 7. Parameters defining the cut wavelengths for the 
order concatenation 



3175 3180 3185 3190 3195 3200 3205 
Wavelength 

Fig. 6. As Figure Q for the LWR camera. In the spectral 
region shown orders do not overlap. The INES concate- 
nated spectrum has no data in that region. 



in the LWR camera, where the S/N ratio in order m is sys- 
tematically higher than in order m-1 in the overlap region. 
In these cases only the points of order m are taken. 

The above defined cut wavelengths (i.e. end wave- 
lengths of order m) can be computed as a function of order 
number as: 



\cut (m) ^ A-\ \ 



(3) 



with the values of A, B and C given in Table 0. 

For non overlapping orders (lower than 73, 77 and 76 
for SWP, LWP and LWR, respectively) only photometri- 
cally corrected pixels have been included in the concate- 
nated spectra (see Fig. ||). The concatenated spectra cover 
the same spectral range as the INES low resolution spec- 
tra, i.e. 1150-1980 A for the SWP camera, and 1850-3350 
A for LWP and LWR. 

Figures ^, [s] and ^ show examples of the concatenation 
procedure for the three cameras. 



Camera 


A 


B 


C 


Orders 


SWP LAP 


24.3952 


132875.4838 


325840.9715 


120-73 


SAP 


22.2095 


133293.4862 


300351.2209 


120-73 


LWP LAP 


-7.9697 


233257.6280 





124-77 


SAP 


-7.7959 


233382.6450 





124-77 


LWR LAP 


-11.3459 


233737.5903 





119-76 


SAP 


-11.2214 


233876.9950 





119-76 



3.1.2. The error vector 

The NE WS IP S piocessmg provides an error vector for the 
high resolution spectra which is computed simply as the 
sum along the extraction slit of the noise values for the 
individual pixels, as derived from the camera noise model. 
Unlike the "sigma" of the low resolution data, the "sigma" 
vector in the MXHI files is not flux calibrated but given 
in FN (Flux Number) units. 

In the INES high resolution data, the "sigma" spec- 
trum is provided in absolute flux units. The calibration is 
performed by applying to the MXHI error vector the high 
resolution calibration and the time sensitivity degradation 
correction. 



3.2. The resampled spectra 

In the INES Archive, each high resolution image has an as- 
sociated "rebinned" spectrum, which is obtained by rebin- 
ning the " concatenated" spectrum at the same wavelength 
step size as low resolution data. This data set represents 
an important complement to the low resolution archive, 
and it is especially useful for time variability studies. The 
rebinned spectra have not been convolved with the low res- 
olution Point Spread Function, and therefore have a better 
spectral resolution than low dispersion spectra. Examples 
of rebinned spectra are shown in Figures m and ra. 

The high resolution concatenated spectra (derived as 
described in the previous Section) are resampled into the 
low resolution wavelength space following the procedure 
detailed below. 



3.2.1. The rebinning procedure 

The concatenated spectra have been resampled into the 
INES low resolution wavelength domain as deflned in Pa- 
per I. The sampling interval is 1.6764 A/pixel and 2.6693 
A/pixel, for the short wavelength and the long wavelength 
ranges, respectively, and the wavelength coverage is 1150- 
1980 A for SWP and 1850-3350 A for LWP and LWR. 
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Fig. 7. Comparison of low resolution (thin lino) and high resolution rebinned (thiek line) spectra of the standard star 
BD+28 4211. The top panel shows the flux spectra. Flagged pixels are marked with diamonds (low resolution) and 
squares (rebinned). The two gaps in the high resolution spectra longward 1900 A correspond to the regions where 
spectral orders do not overlap, which have been assigned zero flux. The inset in the top panel shows in more detail the 
region 1300-1400A. The error spectrum is shown in the bottom panel (thin line: low resolution, thick line: rebinned). 
The characteristic pattern of the errors of the rebinned spectrum is due to the lower signal-to-noise ratio at the edges 
of the individual echelle orders. 



The resampling has been performed so that the total 3.2.2. Errors 



flux is conserved, that is, if n pixels with fluxes /i, /2, 
/„ are rebinned into one, the total flux in the bin is: 



F = ^(A,-A,_i)(/, + /,_i)/2. 

i=l 



(4) 



where the flux at the bin edges {i—1, i~n) is calculated 
by linearly interpolating between the two adjacent pixels. 
The flux of the final pixel is: 



The rebinned error spectrum is computed from the con- 
catenated error spectrum according to the following ex- 
pression: 



E 



VE' 



(6) 



where e^ are the errors of the original pixels in the con- 
catenated spectrum. 



Flux — F / step 



(5) 3.2.3. Flagging 



being "step" the low resolution pixel size defined above. 



The quality fiag assigned to each pixel in the resampled 
spectrum is the sum of the flags of the original high reso- 
lution pixels. Only the most relevant quality flags present 
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Fig. 8. Similar to Figure but comparing two LWP spec- 
tra of Nova Cygni 1992. The errors are shown in the bot- 
tom panel. 

in the high resolution spectrum have been transmitted to 
the rebinned spectrum: 

-8192: Missing minor frames in extracted spectrum 

-1024: Saturated pixel 

-16: Microphonic noise (for the LWR camera only) 

-8: Potential DMU corrupted pixel 

-2: Uncalibrated data point 

Other flags (e.g. reseau marks) have not been taken 
into account to avoid that a too large fraction of the pixels 
in the output spectrum come out flagged with error con- 
ditions, despite their quality is not significantly affected. 
Pixels corresponding to the gaps between non-overlapping 
orders are flagged with "-2" . 

4. The Correction to the Wavelength Scale 

We already pointed out in Section that there is a signif- 
icant discrepancy between the wavelength scales of short 
and long wavelength spectra processed with NEWSIPS. 
This inconsistency, which is well above the accuracy of the 
wavelength calibration, was already present in lUESIPS 
data (Nichols-Bohhn and Fesen 1986, 1990). 



Table || presents a summary of the radial velocities of 
interstellar features in several stars, taken from Table El 
together with values from the literature (v(lit)). The LWP 
and LWR velocities of Table 0, being consistent with each 
other, have been averaged together into v(LW). 

The literature value for RR Tel is from Tackeray 
(1977). The value quoted for ( Oph measured in the opti- 
cal Ca II K and Na I D lines has been taken from Barlow 
et al. (1995). This value is in good agreement with the 
velocities derived from GHRS ultraviolet spectra by Sav- 
age et al. (1992), -14.9 km s"\ and Brandt et al. (1996), 
-15.4 km s^^. The velocities for rj UMa and ( Gas corre- 
spond to measurements of the Ca II K and Na I D opti- 
cal lines reported by Vallerga et al. (1993). The velocity 
quoted for A Lep refers to the optical Ca II doublet (Frisch 
et al. 1990), which presents two components at 2 and 18 
km s~^, which cannot be resolved with lUE. The spectrum 
of HD 93521 presents up to nine interstellar components, 
with the two strongest ones located at, approximately, -10 
and -60 km s"^ (Spitzer and Fitzpatrick 1993). In the lUE 
spectra all these systems appear blended, and therefore, 
as in the case of A Lep, we cannot compare reliably the 
lUE velocities with the optical values. 

According to the data in Table 0, the mean difference 
between long and short wavelength (large aperture) ve- 
locities is 17.7 km s~^, i.e. SWP velocities are systemati- 
cally more negative. The mean difference between the long 
wavelength and the literature values is 8 km s~^. 

A similar test was made to check the consistency be- 
tween the wavelength scale of spectra taken through the 
large and small apertures. Being the number of small aper- 
ture high resolution spectra very limited, useful data were 
available only for the star ( Oph in the SWP and LWR 
cameras. The wavelength scales of the small and large 
aperture spectra are fully consistent in the short wave- 
length range: v(LAP)-v(SAP) = 1.1±6.7 km s^^, while for 
the LWR camera a significant difference is found: v(LAP)- 
v(SAP) = -13.7±4.1 km s~^ The lack of a suitable data 
set precludes an accurate determination of the offset be- 
tween the large and small aperture scales in LWP spectra, 
but the limited tests performed seem to indicate that small 
aperture velocities are systematically lower, although the 
actual difference cannot be quantified. 

The reason for the discrepancy between the short and 
long wavelength range velocity scales is not clear, while the 
large/small aperture discrepancy in LW spectra is most 
likely related to the transfer of the dispersion constants 
from the small to the large aperture: the dispersion rela- 
tions were derived from spectra taken through the small 
aperture and then transferred to the large aperture on the 
basis of the assumed aperture separations. 

In order to provide an internally consistent wavelength 
scale within the INES system, a velocity correction of 
-t-17.7 km s^^ has been applied to the wavelength scale 
of SWP high resolution spectra. The wavelength scale of 
LWP/LWR small aperture spectra has been corrected by 
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Table 8. Mean radial velocities of interstellar lines obtained from NEWSIPS SWP and LWP/LWR large aperture 
spectra 



Target 


v(SWP) 


v(LW) 


v(lit) 


v(LW)-v(SWP) 


v(lit)-v(LW) 




km s~^ 


km s^^ 


km s 


km s^ 


km s^^ 


RRTel 


-69.5 


-50.0 


-61.8 


19.5 


-11.8 


COph 


-24.7 


-11.0 


-14.8 


13.7 


-3.8 


BD-h28 4211 


-22.8 


-3.8 




19.0 




BD-^-75 325 


-16.4 


6.3 




22.7 




HD 60753 


18.8 


31.7 




12.9 




HD 93521 


-38.8 


-20.1 


-10,-60 


18.8 




Tj UMa 




2.4 


-2.8 




-5.2 


A Lep 




20.1 


2,18 






( Gas 




1.2 


-10.3 




-11.5 



Mean difference v(LW)-v(SWP) = 17.7±3.7 km s" 
Mean difference v(Iit)-v(LW)=-8.0±4.2 km s"^ 



+ 13.7 km s ^. With these corrections, the INES velocity 
scale is consistent with the optical determinations. 

5. Conclusions 

In this paper we have discussed the overall quality of lUE 
high resolution spectra processed with the NEWSIPS sys- 
tem. 

The stability of the wavelength scale (« 5 km s^^) 
is within the limits imposed by the acquisition and track- 
ing accuracy. No appreciable distortions in the wavelength 
scale over the full spectral range or during the spacecraft 
lifetime have been found. 

A discrepancy of 9 km s^^ has been found in the ve- 
locities derived from the two components of the Mg II 
doublet at 2800 A in the LWP camera. The correct veloc- 
ity is provided by the K line measured on spectral order 
83. No such discrepancy has been found in LWR spectra. 

The wavelength scales of NEWSIPS short and long 
cameras present an inconsistency, which is well above the 
repeatability errors quoted above. Measurements of nar- 
row interstellar lines have shown that SWP velocities are 
systematically more negative by -17.7 km s~^, on average. 
A similar discrepancy has been detected in long wave- 
length small aperture spectra, whose velocities are more 
negative than those from long wavelength large aperture 
spectra by -13.7 km s^^. 

The determination of the inter-order background has 
greatly improved with respect to the the lUESIPS sys- 
tem, especially at the shortest wavelengths, as shown by 
the absence of negative fluxes in the cores of saturated 
absorption lines and by the greater accuracy of equivalent 
width measurements. 

The INES system derives two spectra from each orig- 
inal NEWSIPS MXHI file. In the first one, the "concate- 
nated" spectrum, the spectral orders are merged together 
and the overlap regions are suppressed according to an al- 
gorithm which computes suitable cut wavelengths which 



maximize the signal-to-noise ratio. This spectrum contains 
an error vector calibrated into absolute flux units, which is 
not available in NEWSIPS data. Since the INES high res- 
olution spectra are obtained from NEWSIPS MXHI files, 
all previous considerations about the stability of the wave- 
length scale, the stability and accuracy of the flux scale 
and the validity of the background extraction are applica- 
ble to them. The second output product is the "rebinned" 
spectrum, which is the "concatenated" spectrum after re- 
sampling at the low resolution wavelength step. 

To correct for the discrepancies found in the NEWSIPS 
high resolution wavelength scale, a velocity correction of 
+17.7 km s^^ for SWP spectra and of +13.5 km s^^ for 
LWP/LWR small aperture spectra has been applied to the 
INES "concatenated" spectra. With this correction, the 
overall INES velocity scale is self-consistent, and agrees 
to within 8 km s^^ with the optical velocity scale. 
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